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ABSTRACT 

Deposition of a massive (10 4 to 10 5 Mq) giant molecular cloud (GMC) into the inner 
parsec of the Galaxy is widely believed to explain the origin of over a hundred unusually 
massive young stars born there ~ 6 Myr ago. An unknown fraction of that gas could 
have been accreted by Sgr A*, the supermassive black hole (SMBH) of the Milky Way. 
It has been recently suggested that two observed 7-ray-emitting bubbles emanating 
from the very center of our Galaxy were inflated by this putative activity of Sgr A* . 
We run a suite of numerical simulations to test whether the observed morphology 
of the bubbles could be due to the collimation of a wide angle outflow from Sgr A* 
by the disc- like Central Molecular Zone (CMZ), a well known massive repository of 
molecular gas in the central ~ 200 pc. We find that an Eddington-limited outburst 
of Sgr A* lasting ~ 1 Myr is required to reproduce the morphology of the Fermi 
bubbles, suggesting that the GMC mass was ~ 10 5 M Q and it was mainly accreted by 
Sgr A* rather than used to make stars. We also find that the outflow from Sgr A* 
enforces strong angular momentum mixing in the CMZ disc, robustly sculpting it into 
a much narrower structure - a ring - perhaps synonymous with the recently reported 
"Herschel ring". In addition, we find that Sgr A* outflow is likely to have induced 
formation of massive star-forming GMCs in the CMZ. In this scenario, the Arches 
and Quintuplet clusters, the two observed young star clusters in the central tens of 
parsecs of the Galaxy, and also GMCs such as Sgr B2, owe their existence to the recent 
Sgr A* activity. 

Key words: galaxies:evolution - galaxies: individual: Milky Way - quasars: general - 
black hole physics - accretion 



_C3_ 1 INTRODUCTION few xlO 39 erg s " 1 in X-rays (e.g. iRevnivtsev et alJliool : 

,„,„ , , iPonti et al.ll2010l ). This variability may reflect f eeding events 

The Galactic centre (GC) is a particularly interesting as- „ „ . 177 — . — 7 — n 

. , f , ? rrom a tew pc-scale molecular gas reservoirs (Morris et al. 

trophysical location, primarily (hie to the presence or the HTT^K • 1 -,-, ■ ,1 ■ 1 r ,1 

f 1999) or variability in the wind capture rate trom the 
closest supermassive black hole (SMBH) - Sgr A , whose . „ . „ , , , _ , . ,, 

, f „, , r-— — tr -ni — ^ 1 ' -1 young massive stars near Sgr A that teed it presently 

mass isMbh ~4x 1O 6 M (|Schodel et aUlgOOa iGhez et all J L, * l l bnosh 



l2005l . l200ah . Although the Soltan relation (|SoltarJll982l )~irn" 



plies that Sgr A* gained most of its mass through luminous A recent observation may shed a rather unexpected 

accretion, it is currently very dim in comparison with ac- light on a much earlier but still rather recent activity of 

tive galactic nuclei (AGN): its X-ray luminosity is less than Sgr A*. Fermi-LAT data has been recently analysed to re- 

~ lCP^Z/Edd (wher e Z/Edd ~ a few x 10 erg s _1 is its Ed- veal two giant 7-ray emitting bubbles, d isposed symme tri- 

dington luminosity; iBaganoff et al.l 120031 ). Sgr A* dimness cally on eit her side of the Gal actic plane (|Su et alj|2010l : al- 

has stimulated development of radiatively inn eficient ac- though see lDobler et al.ll201ll for a different interpretation). 

cretion and accretion/outflow solu tions (e.g. jNaravan et al.l They are roughly teardrop-shaped and extend ~ 8 — 10 kpc 

1 199a : iBlandford &: Begelmanlll999l ). However, X-ray reflec- above and below the plane, but are centred on Sgr A* with 

tion nebulae suggest that Sgr A* might have been much a narrow (d ~ 100 pc) waist along the plane. The limbs 

brighter about a hundred years ago, with luminosity of a of the lobes coincide with the exte nded structure seen i n 

medium-energy X-rays by ROSAT (|Snowden et all ll997J ). 

The lobes have 7-ray lu minosity L-, ~ 4 x 10 37 ergs -1 . Ob- 

* E-mail: kastytis.zubovasaastro.le.ac.uk servational constraints (|Su et al.l 120101 ) and emission mod- 
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elling ([Crocker fc Aharonianll2011f ) allows one to estimate 
the total kinetic and thermal energy content of the Fermi 
bubbles as E buh ~ io 54 " 55 erg. 

There has been a number of suggest ions as t o what 
the origin of the Fermi bubbles might be. ISu et al.l (|2010l ) 
discuss several physical processes and provide a constraint 
that if the bubbles are older than a few x 10 6 yr, the 7- 
ray emission must be powered by ions rather than electrons 
due to a short cooling time of the latter (see their section 
7.1 and Figure 28), unless electrons are continuously acceler- 
ated within the bubbles (which m ay n ot be unreasonable; see 
i]5.4l below). ICrocker et al.l l|201ll ) and lCrocker fc Aharonianl 
(J2011J) detailed these arguments further and suggested that 
the emission is powered by cosmic ray (CR) protons rather 
than electrons. They further consider a quasi-steady state 
model in which the CR protons are continuously injected 
by supernova explosions. CR protons and heavier ions are 
then trapped inside the bubbles for the lifetime of the lat- 
ter (which the authors require to be a bout 10 Gyrs). On 
the other hand, iMertsch fe Sarkar] (|2011T ) argued that the 
emission spectrum of the bubbles is inconsistent with CR 
protons and therefore the bubbles must be a recent, cosmic 
ray electron-powered, feature. 

The scenario with a more recent orig in of the bubbles 
has b een investigated by several authors. iGuo fc Mathews! 
J201ll ') suggest that a jet launched by Sgr A* 1-2 Myr 
ago cou ld create the morph ology and emission structure ob- 
served. ICheng et alj (|2011f ) argue that the bubbles are in- 
flated by episodic Sgr A* activity caus ed by tidal disruptions 
of stars passing too close to Sgr A* . IZubovas et al.l l|201ll . 
hereafter Paper I) argue that an almost spherical outflow 
from Sgr A* caused by a short burst of AGN activity coin- 
cide nt with the well-know n star formation (SF) event 6 Myr 
ago (|Paumard et al.ll2006r ). can be the origin of the bubbles. 
In the model, the outflow from Sgr A* is suggested to be 
quasi-spherical at launch as is required to explain the ob- 
served Mbh — Q" relations for cla ssical bulges and elliptical 
galaxies (|Kmdl2003l . l2005l . l2010ah . The outflow becomes col- 
limated by interaction with a disc-like structure, the Central 
Molecular Zone - a dense ring (aspect ratio H/R ~ 0.2 — 0.3) 
of predominantly molecular gas extending from the Galaxy 
centre to ~ 200 pc. This collimation explains the large 
width of the bubbles and their symmetry with respect to 
the Galactic plane, which is not at all guaranteed by a jet 
model, since jet directions are known to be completely un- 
correlated with the l arge-scale structure of t h e host galaxies 
Kinney et ail |2000| : iNagar fc Wilson! 1 19991 ). IZubovas et all 



20111 ) estimated the duration of the quasar phase to be at 



the minimum t q = 5 x 10 yr and the mass of gas accreted 
by Sgr A*, AM > 4x 1O 3 M , comparable to the mass in the 
young stars in the cen tral parsec, M* ~ few to ten x 10 3 M© 
l|Paumard et alj|2006h . 

In this paper, we test the analytical results of Paper I 
by performing numerical simulations. We broadly follow the 
methodology of simulating spherical out flows from accret- 
ing SMBHs of Navakshin fc Powed |2010l ). These SPH sim- 
ulations involve a prescription for passing momentum and 
energy of the outflow to the ambient gas {Navakshi n et al 



2009). The method has been shown IjNavakshin fc Power 



bles can be reproduced with this model. In addition, we are 
interested in constraints that the simulations could place on 
a key unknown of the problem - the duration of Sgr A* ac- 
tivity. We find that a longer quasar outburst is needed than 
found in Paper I, most likely because the energy produced 
by Sgr A* outflow is spread more widely in the simulations, 
e.g., outside the bubbles themselves, than assumed analyti- 
cally. 

The Paper is structured as follows. In Section [5] wc 
review our analytical model and the main results of Paper I. 
In Section|3j we describe the numerical model and the initial 
conditions of the simulations. We follow this by presenting 
the results in Section [4] and discussing their implications in 
Section [5] We briefly summarize and conclude in Section [6] 



2 A PHYSICAL MODEL FOR THE FERMI 
BUBBLES 

2.1 Sgr A* as the power source 

As explained in Paper I, our Sgr A* feedback model is based 
on the analytical mod els of f ast w i nd o utflows from super- 
massive black holes bv lKind (|2003l . l2005l . l2010bl lah. tailored 
to the present day inner Milky Way. Here we summarise the 
key assumptions behind this model. 

Fast (v ~ 0.1c) wide-angle outflows from AGN have 
been revealed by observations of blue-shifted X-ray absorp- 
tion lines; they a ppear to be a ubiquito us feature, detected in 



ap 



40% of AGN (iTombesi et al|l2010al lbh. Recent XMM ob- 



d. 

n: 



20101 ) to reproduce the analytical r esults of the wind feed- 
back model we use here i|King||2003i ). We aim to investigate 
whether the morphology, size and energy content of the bub- 



servations (|Pounds fc Vaughanl201lf ) of X-ray emission from 
a limb-brightened shell of post-shock gas building up ahead 
of the contact discontinuity corroborate these results and 
the model of King (2003) further. The outflows are probably 
driven by radiation pressure of photons emitted by Sgr A* 
on the gas in the innermost regions of black hole accretion 
flows. While some anisotropy in the outflow may be present 
in general, here we consider a spherical outflow for simplicity 
and demonstrate that it is the properties of the host galaxy 
that may mold the outflows into the teardrop-like shape. 

We assume that Sgr A* accretion rate is mildly super- 
Eddington, and that its luminosity is limited by the Ed- 
dington limit. We do not model the sub-parsec scale disc(s) 
fr om which Sgr A* would be a ccre ting gas here. As show n 
bv lNavakshin fc Cuadral (|2005h and I Alexander et all J2012J ). 
viscous disc time scales are long enough to maintain Sgr A* 
feeding via an accretion disc for at least a million years. Fur- 
thermore, Nayakshin et al (2012, ApJ submitted) show that 
AGN feedback is unable to remove gas discs within the black 
hole's sphere of influence. Physically, this result is easily un- 
derstood: escape velocities for gas in the sub-parsec scale 
disc are much larger than that for gas in the bulge of the 
host, so much more energy needs to be expended to eject 
gas from very near the black hole than from the bulge of 
the Galaxy. This is why a relatively low mass sub-resolution 
reservoir of gas (an accretion disc) could power Sgr A* for 
the duration of the quasar outburst envisioned here, while 
in the meantime Sgr A* drives much more massive outflows 
on scales of the Galaxy's bulge. 

The outflow is believed to be self-regul ating to have 
a scattering optical depth ~ 1 (iKingjuOlOar ). so that the 
photons transfer their momentum to the outflowing wind 
which carries momentum 
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M out v 



(1) 



where v ~ r/c, and 77 ~ 0.1 is the accretion efficiency. We note 
that the mass outflow rate from Sgr A*, M out , is bound to be 
of the order of ~ MEdd, which is negligible on the Galactic 
scales, implying that we can neglect injection of the mass 
into the Galaxy (but not the outflow's momentum or the 
energy, as we discuss now). 



2.2 The effects of the outflow on the diffuse gas in 
the Galaxy 

Since the outflow velocity is very much larger than the es- 
cape velocity from the Galaxy, the outflow propagates with 
an approximately constant velocity until it shocks against 
the ambient gas. For a sufficiently powerful outflow, the 
shocked wind continues to drive an outflow outward, clearing 
the inner regions of the bulge of gas. The Mbh — <? relation 
for elliptical galaxies and classical bulges may have been es- 
tablished by such ou tflows that terminated both host galaxy 
and SMBH growth | |Kinell2005J ). 

The point that we want to emphasize now is that even 
SMBHs lying well below this relation can drive the gas out 
of their host bulges, provided that the gas density is lower 
than t hat expected in th e gas rich era when the bulge itself 
forms (|Kind 120031 . I2005T ) . In our model, gas density in the 
bulge is parametrised through the gas fraction / g : 



2ttGR 2 



2.5 x 10- 26 /-3^kpc g cm" 3 , 



(2) 



where a is the velocity dispersion in the bulge and R is the 
radial coordinate. The oc R~ 2 scaling of the density is nat- 
ural as this is the same as that for the background stellar 
density that dominates the potential in the bulge. We fur- 
ther parametrize / g = 10 _3 /_3, R = _Rk P c kpc and we use 
the value a = 100 km/s, appropriate for th e Milky Way. 
For reference, the model of lKinel (|2003l . I2005T ) assumes that 
f g = 0.16, the cosmological value, during the bulge forma- 
tion epoch. 

The dynamics of the outflow depends crucially on 
whether the shocked wind cools efficiently or not. The shock 
temperature is of the order of 10 11 K, so the only ef- 
fective form of cooling is the inverse Compton effect of 
the shocked outflow's electrons on the Sgr A* radiation 
field (which would dom inate over the whole Galaxy's light; 
ICiotti fc Ostriker]|l997l ). Outside a certain radius, R coo i, the 
radiation field is unable to cool the shock and the out- 
flow transitions from momentum-driven to energy-driven 
(|Zubovas fc King 120121 ) . This radius can be e stimated b y 
comparing the cooling and the flow timescales (|King||2003fl . 
For the parameters of the Milky Way (Mbh = 4 x 10 6 M q ), 
Rcooi is 

R coo i ~ 12 /i'3 2 pc . 



(3) 



Even for the cosmological gas fraction (/_3 = 160), -Rcooi 
is very small compared with the size of the bulge. For the 
problem at hand, where the gas density in the bulge is likely 
to be much smaller than this, this implies that for all the 
relevant parameter space any outflow from Sgr A* quickly 
becomes energy-driven. In this type of outflow, almost all of 
the kinetic power of the outflow, 



-Bou 



-M out v 



^ m - v r 



(4) 



is passed to the ambient gas, heating and pushing it outward. 
The shocked shell rapidly (in t ^ 10 5 yr) attains a constant 
velocity, which for a spherically symmetric outflow is given 
by 

1/3 



2r)cTc 0.16 M BH 

~~ 3 yT^T 



1920a 2 / 3 /r 3 1/3 ' 1/3 



I km s 



.(5) 



where M a is the value of the SMB H mass expected from the 
M -a relation JKing et alj|201ll ). Mbh ^ 0.2M CT and I ~ 1 
is the Eddington ratio. 

Once the quasar switches off, the shell continues to ex- 
pand, but stalls and eventually stops (assuming that the 
isothermal potential of the bulge extends to infinity). The 
maximum radius which the shell reaches is 

install — — Ro — — tq, (6) 

a a 

where t q is the duration of the quasar outburst. 

Using the size of the observed Fermi bubbles 



(R* 



10 kpc) we can constrain the free parameters of 



our model, f-3 and i q . We estimate the gas fraction from 
equation ((3J, requiring that the bubbles must have reached 
the present radius in 6 Myr. This is equivalent to assuming 
that the shell's mean velocity is (v) > 1600 km/s (see Paper 
I for details), and yields 



f- 



/ 1920 \ 

Vieooy 



l~2l 



(7) 



We see that the gas fraction in the Milky Way halo must 
have been similar to 10 -3 in order for the bubbles to be in- 
flated within our model. Numerical simulations below con- 
firm this. 

Requiring that the bubble stalling radius should be 
greater than their current radius, equations © and (JSJ) give 



t q > 2.5 x 10 5 /5 3 3 yr 



(8) 



As the bubbles are likely to be still expanding at the present 
time (see Paper I), this result is a lower limit. We note that 
the estimate given by equation [8] is a factor of 5 larger than 
the estimate obtained in Paper I. As both estimates are 
lower limits on the outburst duration, only the larger one of 
the two is relevant. 

In 0.25 Myr, Sgr A* accreting at its Eddington limit 
consumes AM ~ 2 x 10 4 Mq of gas. Within our feedback 
model, ~2x 10 56 erg is ejected in the outflow's mechani- 
cal energy. This is more than an order of magnitu de larger 
than the bubbles energy content estimated by ISu et al.l 
1)20101) . However, we find that a large fraction of the out- 
flow's energy goes into mechanical work expended to drive 
the ambient cooler mediu m away from the Galaxy's cen- 
tre. |Zubj2vas_&j£ing| (|2012l ) show that even while the quasar 
driving is on, only 1/3 of the energy input is retained in the 
shocked wind while the quasar is active. When the quasar 
turns off, the bubble expands adiabatically, converting its 
thermal energy into kinetic energy of the bubble and the sur- 
rounding shell. Therefore it is likely that the actual amount 
of energy retained by the bubbles is much lower than the 
original energy input by Sgr A* into the outflow, bringing 
the value in line with the observations. Our simulations con- 
firm this prediction; cf. g5.3l 
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2.3 The role of the Central Molecular Zone in 
focusing the outflow 

The shell expansion velocity, v e , depends on the gas density 
in the direction of expansion. Clearly, if the ambient gas dis- 
tribution is not spherically symmetric then the outflow must 
lose its spherical symmetry too. As a minimum effect, the 
velocity of the contact discontinuity must be smaller in the 
directions of denser gas. The most salient feature in the dis- 
tribution of gas in the inner Galaxy is that a good fraction 
of it - mainly the cold molecular gas - lies close to the plane 
of the Galaxy. For what follows it is t he Central Molecu- 
lar Zone (CMZ; lMorris fc Serabvnlll99q ). a massive disc-like 
molecular gas feature in the inner ~ 200 pc, that matters 
the most. As we argued in Paper I, the CMZ presents an 
almost impassable barrier to the Sgr A* outflow. This can 
be seen from the fact that the weight of the CMZ, 

... GM onc McMz 2Mcmzct 34 , . 

Wcmz ^j = — „ ~ 6.5 x 10 dyn, (9) 



its timesteps and move radially with a set constant velocity 
(f w = 0.1c). Each particle carries linear momentum 



r?2 

-"-CMZ 



-RCMZ 



where 7?cmz is the mean radius of the CMZ and M enc is 
the mass enclosed within this radius, is greater than the 
momentum flux impacting on the CMZ from Sgr A* outflow 
(which is the outward force on the CMZ): 

Fcmz ~ = ~ 2 x 10 dyn. (10) 

C itCMZ 

Therefore we expect the CMZ to be a very effective obstacle 
to the outflow (this estimate does not take into account en- 
ergy deposition by the outflow into the CMZ, which makes 
the latter somewhat more prone to the feedback; see H5.6p . 
Physically, we expect Sgr A* outflow impacting the CMZ to 
shock and thermalise. As the CMZ weight is so large, the 
outflow is either completely stopped or nearly stalled in the 
plane of the Galaxy. The shocked outflow gas cannot simply 
pile up there, however. Indeed, if that were the case then the 
pressure (thermal energy density) would increase in that lo- 
cation without limit. The thermal pressure of the shocked 
gas in the directions perpendicular to the Galactic plane is 
much lower because the outflow proceeds in those directions 
easily. Thus there is a strong pressure gradient in the shocked 
outflow gas pointing towards the Galactic plane. This pres- 
sure gradient clearly must launch a "secondary" thermally 
driven outflow away from the Galactic plane, efficiently col- 
limating Sgr A* outflow into these directions. 



3 SIMULATION SETUP 

3.1 Numerical method 

Our workhorse code for solving gas dynamics in the fixed 
potential of the bulge and the black hole is GADGET-3, a n 
updated version of the code presented in ISpringell (|2005h . 
Feedback from the SMBH is implemente d with the 'vir- 
tual p article' method explained in detail in lNavakshin et al.l 
1)20091 ). Rather than hydrodynamically modelling the low- 
density fast wind flowing from the SMBH, its effects on the 
gas distribution at large are simulated via a distinct popu- 
lation of particles, called 'virtual particles' since they exist 
only during their propagation from the source, Sgr A*, to 
the point of interaction with the ambient gas. These par- 
ticles are ejected isotropically from the SMBH at each of 



Pi 



i/EddiBH 

cNl, 



(11) 



and kinetic energy _E 7 = p 7 u w /2, where £bh is the SMBH 
timestep, and divided equally among iV 7 virtual particles 
ejected in that timestep. The scheme is designed to eject 
the correct momentum flux from the SMBH, I/Edd/c, and 
the number of particles ejected in a time step, JV 7 , is selected 
so that the interactions with SPH particles are sufficiently 
frequent to reduce th e (random-number-gen erated) noise to 
acceptable levels (see iNavakshin et al.ll2009l ). 

Each virtual particle has a search radius (a multiple 
of the smoothing length of an SPH particle), within which 
search for potential SPH neighbors is performed "on the 
fly". Whenever an SPH particle kernel contains a virtual 
particle, the two are considered interacting, and a fraction 
of the virtual particle's momentum and energy is transferred 
to the gas. When the virtual particle energy drops below 1% 
of its original energy, the particle is destroyed. 

The kinetic energy of the virtual particles is passed to 
the SPH particles if the shocked wind does not cool, i.e. 
if R > i?cooi (eq. |3]) . We have implemented a step function 
transition between the pure momentum feedback regime and 
the energy feedback one for this. However, testing showed 
that since i? CO oi is very small in the case considered here, 
there is no significant difference in results between simula- 
tions utilising the switch and those employing the energy 
feedback scheme throughout the simulation, so we use the 
latter below. 

While we do not model hydrodynamically the reverse 
wind shock, which thermalizes the wind and creates a pres- 
surized bubble that can expand adiabatically, we find that 
this process is somewhat mimicked in the simulations by low 
density gas that is present in the otherwise evacuated bub- 
ble. These particles predominantly originate on the surfaces 
of the Central Molecular Zone (see Sections l3.2l and l4.1.3l be- 
low), where they are heated to very high temperatures and 
rise to fill the voids; subsequently quasar wind heats these 
particles even further. The bubbles are found to be signifi- 
cantly overpressurized with respect to the surrounding gas. 
This allows the bubble to expand thermally as expected if 
the reverse shock were modelled by SPH (see Section 14.11 
below). 

In all the simulations presented below, we use an ideal 
equation of state for the gas. The gas pressure is given by 
P — pkT I [i, where mean molecular weight, p = 0.63m p (as- 
suming ionised gas of Solar abundance), k is the Boltzmann's 
constant, and p and T are the gas density and temperature, 
respectively. In all simulations except for one, we employ 
the standard GADGET optically t hin c ooling prescription 
based on the lSutherland fc Dopital (|l9931 ) cooling curves, in 
addition to heating from the virtual particles on the contact 
discontinuity of the bubble. In the only exception to this, 
simulation 'Cool' (see Table [T}, we check the sensitivity of 
our results to the assumed cooling function by utilising the 
optically thin radiative cooling rates for gas ionised (and 
also heated) by the quasar radiation field as calculated by 
iSazonov et al.l (2005), who also provided an analytical fit to 
the respective rates that we use here. 
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Test 


u 


-* *p,halo 


tq/Myr 


\ R J cmz 


Afcm Z /M0 


iVp^mz 


Rb/kpc 


db/kpc 


h b /kpc 


«v/km s 1 


tJh/km s 1 


Base 


10~ 3 


7.5 x 10 5 


1 


0.25 


10 8 


10 6 


11.5 


9 


~ 1 


1090 


670 


HR-low 


10~ 3 


7.5 x 10 5 


1 


0.125 


10 8 


10 6 


11.5 


9 


~ 1 


1110 


680 


Cool* 


10~ 3 


7.5 x 10 5 


1 


0.25 


10 8 


10 6 


9.5 


10 


~ 1 


820 


650 


Fg-low 


4 x 10~ 4 


3.0 x 10 5 


1 


0.25 


10 8 


10 6 


~ 15 


~ 10 


2.5 


2560 


610 


Fg-high 


4 x 10~ 3 


3.0 x 10 6 


1 


0.25 


10 8 


10 6 


6 


6 


~ 1 


220 


100 


Tq-low 


10~ 3 


7.5 x 10 5 


0.3 


0.25 


10 8 


10 6 


7 


4 


2 


300 


360 


Both-low 


4 x 10~ 4 


3.0 x 10 5 


0.3 


0.25 


10 8 


10" 


~ 7 


~ 4 


1 


960 


270 


Mc-low 


10~ 3 


7.5 x 10 5 


1 


0.25 


10 7 


10 5 


7 


8 


2 


330 


170 



Table 1. Simulation parameters and main results. From left to right, the parameters are: Test ID, gas fraction, number of SPH particles 
in the halo, quasar outburst duration, CMZ mass, CMZ scale height, number of particles in the CMZ. The results are: bubble height, 
width and distance between its lower edge and the SMBH; velocity of the swept-up ISM in the z direction at x = y = and velocity of 
the swept-up ISM in the xy plane at the mid-height of th e bubble; all five at 6 Myr. 
* - Simulation 'Cool' is identical to 'Base', but includes a lSazonov et al.l q2005i ) heating-cooling prescription. 



Simulation snapshots showing density and temperature 
ar e plotted using an an g le slic e projection method presented 
in iNavakshin fc Power! (|2010T ). Specifically, the gas column 
density projected over the y coordinate is calculated by 

£(*,*)= f p(x,y,z)dy, (12) 

where the limits of the integration are given by y(x, z) — 
rtan£ and r = \/x 2 + z 2 . The angle £ is chosen so that 
tan£ = 1/4 throughout this paper. This projection method 
conveniently allows us to get an unobscured look into the 
inner parts of the simulation and yet have enough particles 
at the outer edges for a statistically meaningful figure to be 
derived when plotting. In addition, since most of the results 
presented in this paper are symmetrical around the z axis, 
the snapshot plots, where applicable, are divided vertically 
and show the surface density in blue-white on the left and 
temperature in red-orange on the right. 

3.2 Galaxy model and initial conditions 

The initial setup for the Galaxy consists of three compo- 
nents, described below and summarized in Table [T] The 
'Base' model is the model which be st reproduces the Fermi 
bubble observations (|Su et alj|2010h . We analyze this simu- 
lation in detail below, and study the robustness of our results 
by varying the free parameters of the model. 

The whole computational domain is embedded in a 
static isothermal background potential with a — 100 km/s. 
In the centre of the coordinate system, fixed in space, is the 
SMBH. While we do not model accretion onto Sgr A* in de- 
tail (cf. H2.1|l . ocasionally SPH particles may get very close 
to the SMBH, especially closer to the end of the simulations 
when feedback from Sgr A* is turned off. To avoid very small 
time steps, and thus very high numerical costs, associated 
with these "uninteresting" particles, we remove them using 
the sink particle formalism if the SPH particles are closer 
than r s i n k = 0.1 pc to Sgr A*. As explained above, Sgr A* 
is "turned on" at the start of each simulation, radiating at 
its Eddington limit for a duration t q , which is a free param- 
eter of the model. Most of the models presented here use 
t q = 1 Myr (since this was found to give the best results), 
although we explore shorter outbursts as well. 



As mentioned above, we position a mas sive disc of 
gas - the Central Molecular Zone (CMZ; lMorris fc Serabvnl 



Il996l ) - in the plane of the Galaxy at a distance between 
ftn,cmz = 5 pc and .Rout, cmz = 200 pc from Sgr A*. The inner 
radius is chosen so that the whole CMZ would be comfort- 
ably outside the sphere of influence of Sgr A*(7?j n fl ~ 2 pc), 
but still account for the yet smaller circumnuclear disc of 
gas, which i s observed to lie between ~ 2 and ~ 10 pc from 
the centre l|Guesten et all Il987l : iMorris fc SerabvrJ Il996h . 
Our model CMZ has uniform temperature throughout, cho- 
sen to give the appropriate constant scale height aspect ratio 
(H/R — 0.25 or 0.125, with a corresponding difference in the 
adopted temperature floor for the simulations). It is rota- 
tionally supported in the radial direction, with v TOt — y/2a. 
The radial density distribution follows a p oc R~ 2 power law, 
and the mass of the disc is set to M cmz = 1O 8 M0 (approx- 
imately the upper limit from observational constraints), al- 
though we ran one simulation with a much lower CMZ mass. 
Despite the high mas s of the CMZ, i t is marginally stable to 
self-gravity, since the lToomrel (1 19641 ) parameter Q ~ 3.4 for 
H/R = 0.25, and 1.7 even for H/R = 0.125. We note that 
we do not fine-tune the CMZ to be only marginally sta- 
ble in the simulations; it is a natural outcome of using the 
observationally constrained CMZ parameters. This result is 
probably not a coincidence: massive cold dis cs are widely be- 
lieved to be self-regulating to have Q ~ 1 (|Goodmanl 120031 ; 
iThompson et al.ll2005l ; INavakshin et al.ll2012i ). In principle, 
such a low Q parameter might lead to fragmentation. How- 
ever, since we set the temperature floor in the simulations to 
be equal to the initial CMZ tempe rature, we do not expect a 
significant fragmentation to occur !|Gammiell200ll ; lRice et al.l 
120051 ). This assumption is consistent with simulation results 
(see Sections 14. 1.3l and 15. 6l below~). The SPH particle mass of 
the simulations is M part = lOOM© , giving a minimum mass 
resolution of ~ 4OOOM . 

Finally, there is a diffuse spherically symmetric gaseous 
"halo" extending between r ln = 1 pc and r out = 15 kpc with 
density following the isothermal profile given by equation [5] 
The gas fraction / g is a free parameter of our model and is 
varied between 4 x 10 -4 and 4 x 10~ 3 in different tests, with 
the fiducial value of 10 -3 for the 'Base' model (this corre- 
sponds to a particle density 6 x 10 -4 cm -3 at R = 5 kpc; see 
the Discussion section for the validity of this choice for gas 
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Figure 1. Left: The face-on column density of the CMZ disc at time t = 0.1 Myr for the 'Base' simulation. Note that the innermost 
region has been partially evacuated by Sgr A* feedback. Right: Cross-sectional plot of gas surface density (left half of the panel, blue- 
white) and temperature (right, red-orange) for the same snapshot. The CMZ (cold dense wedge in the Galactic plane) strongly collimatcs 
the outflow; even though its surfaces are ablated, teardrop-shaped cavities form readily. 



density) . The number of particles in the halo is set by their 
mass (we use same SPH particle masses for the CMZ disc 
and the halo). The gas temperature in the bulge halo is ini- 
tially set to T ha i = T vir = 2.5 x 10 5 K, which corresponds to 
the virial temperature of the bulge. The initial gas velocity 
in the bulge is set to zero. 

This halo setup we use is probably oversimplified. For 
example, we do not account for the likely anisotropy of the 
initial gas distribution of the "halo" due to overall rotation 
of the Galaxy. We address this point qualitatively in the 
Discussion section. 



4 RESULTS 

Table Q] shows the list, parameters and main results of all of 
the simulations that we present in this paper. Results of the 
three simulations in the top of the Table, separated from the 
rest by a horizontal line, appear to b e a reasonable match 
to the Fermi Bubble observations bv lSu et al.l (|2010r) . The 
rest of the simulations produce bubbles that are unlike the 
observed ones. 

We first describe the evolution and properties of the 
'Base' model (/ g = 10~ 3 ,t q = 1 Myr). 



4.1 Base simulation 

4-1.1 Small and intermediate scales 

Overall, the dynamics of gas in the simulation closely follows 
our analytical expectations. As soon as the quasar switches 
on, the spherically symmetric outflow hits both the CMZ 
and the halo gas. The left panel of Figure Q] shows the face-on 
view of the CMZ disc at time t = 0.1 Myr. Note that only the 
central ~ 25 pc of the disc were evacuated by outflow from 
Sgr A*. The right panel of the same figure shows the edge- 
on projections of both gas column density and temperature 
(which we present in a single panel because of the azimuthal 
symmetry of the gas flow). The arrows show gas velocity 
vectors projected on the plane of the figure. 

As expected, the spherical "halo" (the diffuse gas com- 
ponent) is affected by the outflow much more than the CMZ, 
with the contact discontinuity between the wind and the 
shocked ambient medium at a distance of about 300 pc. We 
also note that the bubbles do contain some gas mainly closer 
to the interface with the upper layers of the CMZ. That gas 
is heated to temperatures above 10 K. There is a transi- 
tion region between the almost-spherical outflow perpendic- 
ular to the Galactic plane and the stalled outflow against 
the CMZ, completing the figure-8 morphology of the whole 
flow. 

As time proceeds the inner hole in the CMZ disc grows 
in size, engulfing most of the disc by the end of the simula- 
tion (see 34.1.31 for a fuller discussion of this) . The cavities 
opened by the outflow in the directions perpendicular to the 
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Figure 2. Left: Side view of gas surface density (left half of the panel, blue-white) and temperature (right, red-orange) in the central 
2 kpc of the 'Base' simulation at t = 1 Myr. The CMZ has been enveloped by the external ISM shock fronts, but it still maintains 
the strong collimation of the diffuse cavities. The cavities are also filled with hot 10 8 — 10 9 K gas, which allows them to expand in all 
directions once the feedback from Sgr A* has switched off. Right: Central 500 pc of the 'Base' simulation at t = 6 Myr. The CMZ 
remains, although its density structure is perturbed. There is also a "back flow" of warm (T ~ T v ; r ) gas into the central regions, evacuated 
by the buoyant rise of the bubbles. 



Galactic plane grow even more. The left panel of Figure [2] 
shows the edge-on view of the inner 2 kpc of the Galaxy 
at time t — 1 Myr. One notes the hour-glass shape of the 
cavities opened by the outflow, and a strong gradient in the 
outflow velocity with angle 8 (tanr? = |z/a:|) meausured from 
the z-axis. In particular, the maximum velocity is reached 
at 6 — 0, such that velocity there is consistent with the ana- 
lytical prediction of v c ~ 2000 km/s, and the minimum is at 
9 = 90° . There is also a "failed" outflow around 6 « 65° ; gas 
flowing along these directions eventually falls on the Galactic 
plane, shadowed by the CMZ from further Sgr A* feedback. 
Some material is ablated from the surfaces of the CMZ. It 
expands and contributes to the tenuous gas filling the cavi- 
ties. 

The right panel of Figure [2] presents the edge-on view 
of the central 500 pc at the end of the simulation, at t = 
6 Myr. Due to the relative buoyancy of the hot gas in the 
cavities with respect to the cooler "ambient" gas, the former 
leaves the region by that time, being replaced by the latter. 
This "back flow" of warm (T ~ T v i r ) gas returning to the 
central region after Sgr A* switched off is clearly seen in the 
pattern of velocity vectors on the left side of the right panel 
of Figure [2] We also point out that while the CMZ has been 
significantly affected by the outflow from Sgr A* , most of it 
remains in the region in the form of a ring discussed further 
in 34X31 



4-1.2 Large scales 

We now discuss the larger scales of the "Base" simulation. 
Figure [3] shows the edge-on views of the simulation domain 
at times t — 1, 3 and 6 Myr (left, middle and right panels, 
respectively). The left panel in particular shows that by the 
time Sgr A* switches off the cavities are still rather small on 
Galactic scales: their height is R ~ 3 kpc and the maximum 
width is d ~ 2.5 kpc. 

Further evolution of the hot bubbles is driven by the in- 
ertia of the outflow, the buoyancy of the bubbles and the fact 
that they are significantly over-pressurized with respect to 
the ambient medium (cf. the left panel of Figure 2]) . Bubble 
expansion proceeds in an almost self-similar fashion, except 
for a ripple at roughly the middle of the bubble height (see 
Fig. [3 middle panel). It is most likely a Kelvin-Helmholtz 
unstable mode which arises due to the material inside the 
bubble moving parallel to the surface of the contact discon- 
tinuity. The ripple rolls over and disappears by t — 6 Myr 
(Fig. S right). 

By the end of the simulation, the cavities have expanded 
to reach a height of ~ 11.5 kpc (Fig. [3] right panel), consis- 
tent with the o bserved vertical extent of the Fermi bubbles 
(|Su et al.ll20lbT ). The width, at ~ 9 kpc, is slightly larger 
than observed (d b s ~ 6 kpc); we return to this point in 
the Discussion section. The cavities are filled with very hot 
(T bub ~ 2 x 10 8 K ~ 17 keV) and diffuse (/ g ,bub ~2x 10" 5 ; 
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Figure 3. Gas evolution on large scales in the 'Base' simulation; the three panels correspond to t = 1,3 and 6 Myr, from left to right. 
Only the positive-z side of the computational domain is shown, due to symmetry around the Galactic plane. The CMZ strongly collimatcs 
the outflow and allows the formation of a teardrop-shaped cavity with a morphology very similar to that of the observed Fermi bubbles. 
The bubbles continue to expand and rise due to high pressure and low density, even once the feedback has switched off. 
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Figure 4. Left: Gas pressure as a function of cylindrical radius (r 2 j = x 2 + y 2 ) for a selection of SPH particles at 1 < \z\ < 2.5 kpc 
(middle of the bubble) for the 'Base' simulation at t = 1 Myr; the solid line is the mean value at each radius. The central cavity is 
significantly overprcssurised in comparison to the ambient medium, leading to lateral expansion of the bubbles. Middle: Same but at 
time t = 6 Myr when the bubble is much larger. The vertical cut chosen is now 4.5 < \z\ < 6 kpc. Note that the gas pressure has dropped 
significantly and is slowly varying accross the bubble and the shocked region. Right: Temperature (green curve and dark blue points, 
left scale) and gas fraction (red curve and light blue points, right scale) against cylindrical radius for the same time and particle cut 
as the middle panel. The hot diffuse inner cavity and dense surrounding medium are obvious and have rather sharp edges (thickness 
~ 1 kpc). Solid lines show mean values of gas fraction and temperature at each radius. 



see Fig. |H right) gas. For completeness, Figure [5] shows the 
time evolution of the bubble's height, R, and width, d. The 
lateral expansion is initially somewhat faster than the ver- 
tical, possibly because some of the material expanding out- 
ward along the disc plane was blown off the top of the CMZ 
disc (rotating gas is easier to blow away due to centrifugal 
force). As is seen in the figure, in the next few Myr, the 
bubbles are as wide as they are tall at t ~ 2 Myr. Subse- 
quently, vertical expansion becomes faster than lateral one. 
This may be driven by buoyancy: as the bubbles continue to 
rise up, their lower edges "lift off" from the plane, and the 
cooler gas can start flowing back along the plane. 

The bubbles are slightly detached from the very centre 
of the Galaxy; the gradual increase in this detachment is 



visible in Figure [3J and also in the right panel of figure [5] 
This detachment is caused by the bubbles rising buoyantly 
out of the Galaxy potential. By t — 6 Myr, the gap between 
the centre and the lower edges of the bubbles is ft ~ 1 kpc, 
in agreement with obser vations that pe rmit any value of the 
gap size below ~ 2 kpc (|Su et al.ll201fj ). 



The left panel of Figure [4] shows the SPH particle pres- 
sure versus their cylindrical radius defined as r 2 yl = x 2 + y 2 
selected at a slab of gas with 1 < \z\ < 2.5 kpc at t — 1 Myr, 
roughly corresponding to the midplane (i.e. half-height) of 
the bubbles at that age. We observe that the difference in 
pressure within the bubble and outside is greater than two 
orders of magnitude. This excess pressure continues to drive 
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Figure 5. Height (solid line) and width (dashed line) of the bub- 
bles as function of time in the 'Base' simulation. 



the expansion in directions both parallel and perpendicular 
to the Galactic plane long after the quasar has switched off. 
The middle panel of figure [4] presents the same as the 
left figure but at the end of the simulation and for a slab 
of material at a larger \z\ (4.5 < \z\ < 6 kpc, corresponding 
roughly to the middle of the bubble that has now risen fur- 
ther from the Galactic plane at t — 6 Myr). Note that by 
t = 6 Myr, the gas pressure inside the bubble has dropped 
by two orders of magnitude; the pressure of the surrounding 
ISM is also lower, but only by about a factor 10, so the pres- 
sure difference between the bubble and its surroundings is 
much smaller. Consequently, the bubble expansion is much 
slower, as seen in Fig. [5] and the right panel of Fig. [3] The 
vertical expansion of the bubbles persists for longer due to 
residual momentum in the gas and buoyancy (see the last 
two columns in Table [TJ. 



4.1.3 Feedback effects on the CMZ 

In the simulation, we find that the CMZ is not dispersed by 
the outflow, a result consistent with the prediction of the 
analytical argument (cf. Section [2.3[) . However, the quasar 
wind is powerful enough to displace the inner ~ 120 pc of 
the CMZ to larger radii, resulting in formation of a dense 
thin ring (Figure [6j left). 

The average radial expansion velocity of the inner parts 
of the CMZ, u r ,cmz ~ 100 km/s, is larger than the sound 
speed in its gas (c s ~ 20 — 40 km/s, for H/R = 0.125 and 
0.25 respectively) and comparable to its rotational velocity. 
Therefore some of the CMZ gas is shock-heated to T ~ 
10 6 K and expands vertically. This allows the quasar wind 
to ablate the outer surfaces of the CMZ further. Gas from 



these regions fills the voids in the halo (cf. Sections l4.1.1l and 
14.1.21 above). This process, however, removes only a small 
amount of mass: by the end of the simulation, the CMZ 
mass (defined as the gas mass within a radius 250 pc in the 
Galactic plane and within z ± 100 pc) has decreased by less 
than 5%. 

We would like to note another more significant dynam- 
ical effect of the Sgr A* outflow on the CMZ that re-shapes 
the initial disc configuration into a ring-like one. The out- 



ward force and pressure of quasar wind cause a significant 
radial mixing in the CMZ. For example, some of the CMZ 
gas on the inner face of the disc is sent outward by the wind 
just above the surface of the outer regions of the disc, but 
stalls later on and falls back on the CMZ in the outer re- 
gions. This failed outflow has a small angular momentum 
in comparison to gas in the outer disc. The angular momen- 
tum of the CMZ gas is then well mixed up. This establishes a 
narrow distribution of specific angular momentum - a ring. 
We note that this mecha nism of ring f orma tion is similar 
in spirit to that found bv lrlobbs et alj (|201ll ) in their sim- 
ulation S30, although there shocks between material with 
different angular momentum were due to initial conditions 
in the collapsing gas shell rather than quasar feedback. 

Subsequently, the ring slowly relaxes and spreads back 
somewhat into a disc configuration due to viscous stresses; 
however, the viscous timescale for this to happen at R = 
100 pc is t viB C ~ 10 8 (0.1 /a) yr for H/R = 0.25, where a is 
the standard IShakura fc Sunvaevl (119731 ) viscosity parame- 
ter. This is much longer than the timescales we are interested 
in, and thus unsurprisingly the ring persists to the end of 
the simulation (Figure [6] right). 

Since the disc is only marginally gravitationally sta- 
ble at the start of the simulation, it comes as no real sur- 
prise that the ring becomes more unstable than the origi- 
nal disc. As a result, some spiral density waves are visible 
in the middle panel of Fig. [5] Furthermore, a very mas- 
sive dense clump forms at t ~ 3.5 Myr and there are hints 
of another forming by the end of the simulation at a po- 
sition {x,y} ~ {80,-20} pc. The radius of the clump is 
r c i ~ 5 pc and its mass m c \ ~ 10 7 Mq. Its density is then 
n c \ ~ 7x 10 5 cm -3 , ~ 500 times greater than the background 
potential density at the clump's radial distance R ~ 100 pc 
and similar to that of dense star forming molecular cloud 
cores. Due to the adopted temperature floor and lack of res- 
olution in the simulations, the cloud cannot fragment into 
smaller globules, but we obviously expect such an object to 
be unstable to gravitational collapse. The result would pre- 
sumably be a massive star cluster. The orbit of the clump 
around Sgr A* is mildly eccentric (e ~ 0.2). This may be 
interesting as a potential route to formation o f the Arches 
clust er which has a rather non-circular orbit (jStolte et al.l 
I2OO81 ) in particular, but it may be also relevant to the origin 
of other young Galactic Centre star clusters and GMCs (cf. 
further discussion in ^5.6.3|) . 



4.2 Dependence on the ambient gas density 

The ambient gas density before Sgr A* quasar outburst is a 
free parameter of the model. Therefore, we varied f g to see 
how our conclusions depend on this paerameter. A fourfold 
increase in f s (Simulation 'Fg-high', Fig. [7] left) results in 
a reduction of the bubble height by ~ 40% for the same 
£ q . This is consistent with the analytical prediction from 
equation © and can be understood in terms of the larger gas 
weight, which for an isothermal distribution is independent 
of radius and is 



W R = 



4/ g cr 4 



(13) 



(King 20031. l2010aT ). Denser gas weighs more, therefore re- 
quiring more energy to be lifted to the same height. Since 
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Figure 6. Density of gas in the central 250 pc of the 'Base' simulation at t = 1, 3 and 6 Myr (left, middle and right panels, respectively), 
as seen from above the Galactic plane. The inner regions of the CMZ are strongly affected by the quasar activity, forming a ring (left 
panel) which persists for several Myr (middle panel) while slowly spreading radially and developing unstable spiral filaments (right panel). 
A self-gravitating clump forms at t = 3.5 Myr, with another about to form by the end of the simulation. The total mass of the CMZ 
remains almost constant throughout the simulation. 



the energy input is the same in both cases, higher values 
of / g result in lower height reached by the bubbles. In fact, 
in simulation 'Fg-high', formally, the stalling radius of the 
bubble is -Rstaii ~ 10 kpc, similar to the size of the observed 
bubbles. However, the stalling time is tstall ~ 50 Myr, which 
is much greater than the time since the hypothesized Sgr A* 
outburst. This shows that relatively high values of / g are def- 
initely disfavored within our model of a recent outburst ori- 
gin for the Fermi bubbles (unless perhaps £ q is much longer 
than a 1 Myr; but see Section [4. 3p . 

Conversely, a simulation with a lower gas density ini- 
tial condition ('Fg-low', Fig. [7] middle) produces a bubble 
that is ~ 40% taller, but its interior temperature is rather 
low, blurring the distinction between the bubble and its 
surroundings. Furthermore, the lower density bubble rises 
further from the Galactic centre, producin g a detachmen t 
hb > 2 kpc, inconsistent with the data (|Su et al.l |2010|) . 
These disagreements allow us to exclude the possibility of a 
significantly lower gas fraction in the halo as well. 

Our simulations show that the width of the bubbles, 
on the other hand, is almost independent of the ambient 
gas density. This is probably because the lateral expansion 
of the bubble is governed not only by the thrust from the 
outflow, but also by the pressure balance on both sides of 
the bubble (Fig. [3] left and middle panels). The ambient 
gas pressure increases with higher / g , but so does the pres- 
sure inside the cavity due to lower bubble volume. The two 
changes compensate for each other, and the pressure in- 
side the bubbles remains several orders of magnitude greater 
than in the external medium, leading to lateral expansion. 
Additionally, this expansion happens on an approximately 
dynamical timescale, which is given by the properties of the 
background potential and is also independent of / g . 



4.3 Dependence on the outburst duration 

We now explore how simulation results depend on the out- 
burst duration. Simulation 'Tq-low' (see Table 1) is identical 
to the 'Base' simulation except £ q = 0.3 Myr. As expected, 
lower tq results in physically smaller bubbles (see Fig. [7] 
right panel). Due to a smaller amount of energy injected 
and a longer duration of adiabatic expansion phase after 
quasar switch off, the contrast between the temperature and 
density in the bubble interior and the surrounding shell is 
also smaller. The bubbles do have lower densities (effective 
/ g ~ 4 x 10~ 4 , compared with the ambient / g ~ 10~ 3 ) and 
somewhat higher temperatures (Tbub ~ 10 7 K) than the sur- 
roundings, and thus do rise buoyantly away from the centre, 
resulting in large (> 3 kpc by 6 Myr) separations of the 
bubbles from the centre. Vertical expansion of the bubble, 
however, is slower than in the 'Base' simulation; the outer 
edges of the bubbles reach only ~ 5 — 6 kpc by the end of 
the simulation. In addition, the bubble width is significantly 
smaller than in the 'Base' simulation, since a lower energy 
input results in a lower pressure inside the bubble, therefore 
lateral expansion is also slower. 

This result provides a rough lower limit to the quasar 
outburst duration required to inflate the bubbles: i q > 
0.3 Myr, and quite likely i q ~ 1 Myr. This is compatible 
with the analytical results, where we found i q > 2.5 x 10 yr 
from morphological arguments. We discuss this point further 
in the Discussion section. 

Given that a lower value of t q produces smaller bubbles, 
while a lower value of / g increases their size (cf. Section ^. 2[) . 
is it possible that reducing both t q and / g may yield as good 
or better morphological fit to the observed bubbles as the 
"Base" simulation? We run a simulation to test this idea (see 
simulation "Both-low" in Table [TJ , but find results incom- 
patible with observations. Figure|9]shows the surface density 
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Figure 7. Density and temperature plots at t = 6 Myr of simulations with varying parameters. Left: higher gas fraction / g = 4 X 10 
('Fg-high'). Middle: lower gas fraction / g = 4 X 10 -4 ('Fg-low'). Right: lower quasar outburst duration t q = 0.3 Myr ('Tq-low'). 
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Figure 8. Density and temperature plots at t = 6 Myr of simulations with varying parameters. Left: CMZ mass decreased by a factor 
10 ('Mc-low'). Middle: CMZ aspect ratio decreased to H/R = 1/8 ('Hr-low'). Right: a different heating-cooling prescription has been 
used ("Cool'). 



and temperature projections at the end of the simulation. 
While the height of the bubbles is comparable to that found 
in the "Base" simulation, the width is too narrow. Also the 
gap between the Galactic centre and the lower edge of the 
bubbles is too large, and the morphology of the whole region 
appears too non -uniform compared with the observations by 
ISu et all (|201Ch . 

Increasing both t q and / g simultaneously by a factor 
larger than a few is also not a viable option. Morphologi- 
cally, had the CMZ been able to withstand onslaught from 
Sgr A* feedback for longer, the dynamics of the gas may 
have produced reasonably shaped hot bubbles. But we are 
already using a CMZ mass that corresponds to the upper 



limits derived from observations IjMorris fc Serahvnl 119961 ). 
and even with £ q = 1 Myr the CMZ is significantly affected 
in the "Base" simulation. It appears to us that increasing £ q 
further may do too large a damage to the CMZ, i.e., drive 
most of it to much larger distances from the Galactic Centre 
than it currently observed. Therefore we feel that tq longer 
than ~ 1 Myr are not very likely, although there remains 
a possibility that the current CMZ is only a remnant of an 
even more massive gas disc that preceeded the Sgr A* quasar 
outburst. 
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Figure 9. The structure of the bubbles at the end of the sim- 
ulation 'Both-low' in which both the duration and background 
density of gas are reduced from that of the 'Base' simulation to 
t q = 0.3 Myr and / g = 4 X 10~ 4 . The morphology of the bubble 
produced in this simulation is inconsistent with the observations. 



4.4 Dependence on the CMZ properties 

Having checked the effects of varying the halo and quasar 
parameters, we now vary the properties of the CMZ. In 
simulation 'Mc-low', we reduce the CMZ mass by a fac- 
tor of 10, making it smaller than the current observational 
estimates (M cmz ~ 3 - 5 x 1O 7 M ; iDahmen et~al] ll998J ; 
iPierce-Price et al.ll2000T ). 

The results are shown in Fig. [HI left. Comparing the fig- 
ure with the right panel of Figure [3j we note that the shape 
of the bubbles and the shocked region is closer to spherical 
in the case of a less massive CMZ. Although the lower CMZ 
mass still yields a weight ~ 5 times larger than what the 
quasar outburst should be able to lift, the vertical density 
stratification in a homogeneous disc makes a large part of 
the CMZ diffuse enough to be blown away. The rest of the 
CMZ is shock-heated to higher temperatures than in the 
'Base' simulation and begins to expand vertically, provid- 
ing more material that can be removed by the wind. As a 
result, most of the CMZ is blown away by ~ 0.7 Myr in sim- 
ulation 'Mc-low', and the outflow proceeds quasi-spherically 
for the last 0.3 Myr of the quasar activity. Although almost 
all of this material eventually accretes back onto the reform- 
ing CMZ disc after the quasar has switched off, the effects 
on the diffuse "halo" gas are profound. The bubble height 
is ~ 40% lower, their width is ~ 20% smaller and the bub- 
bles are significantly detached from the Galactic centre. This 



simulation does not appear to match the morphology of the 
Fermi bubbles as well as the "Base" simulation does. 

We next check the influence of the CMZ geometry by 
reducing its scale height to yield the geometrical aspect ratio 
of H/R — 0.125, with a corresponding decrease in the tem- 
perature floor of the simulation ('HR-low', Fig. [5J middle). 
This value is smaller than the CMZ aspe ct ratio favored by 
the current observation s (H/R ~ 0.15; IPierce-Price et al.l 
120001 ; Ijones et alJl201lT ). so that we effectively cover all the 
reasonable parameter space in terms of H/R for the CMZ. 

We find no qualitative difference between simulations 
'HR-low' and 'Base'. This suggests that even if the exact 
geometry of the CMZ plays a role in determining the bubble 
shape, the magnitude of this effect is quite limited. In all 
simulations, the bubbles are collimated much more strongly 
than pure shielding by the CMZ would suggest (the opening 
angle of the bubbles is fib ~ 0.4 x 47r, while the solid angle 
not obscured by the CMZ in the 'Base' simulation is fi cm z ~ 
0.8 x 4tt), therefore this lack of difference is not particularly 
surprising. 

4.5 Effect of the heating-cooling prescription 

Simulation 'Cool' is identical to the 'Base' simulation, ex- 
cept that we use a physically motivated optically th in quasar 
heating-cooling prescription (ISazonov et al.l 120051 ) instead 
of the standard ISutherland fc Dopital ( 1993T ) one, as in the 
other simulations in this paper. The prescription is based 
on a fit to the radiative heating and cooling rates of the 
gas illuminated by a typical quasar radiation field. The ra- 
diative processes include photoionization heating, Comp- 
ton and inverse-Compton processes, bremsstrahlung and line 
cooling, and assume an optically thin medium, which is well 
justified for the low column densities that we find in this 
paper, E < 10~ 4 g cm -2 . The only region where gas may 
become optically thick to X-rays from Sgr A* is the mid- 
plane of the CMZ disc. However, we find that the higher 
density gas is able to cool efficiently anyway, thus staying 
close to the imposed temperature floor of T = a few times 
10 4 K in any event. 

The overall effect of changing the heating-cooling pre- 
scription on the large scale gas distribution is small (Fig. [5] 
right). The bubble interior is ~ 50% hotter (T bub ~ 3x 10 8 K 
as opposed to 2 x 10 8 K in the 'Base' simulation), most likely 
due to less efficient cooling in the adopted prescription at 
t > t q . The temperature of the outer shell is practically 
the same. Morphologically, this difference results in bub- 
bles that have ~ 10% lower height and practically the same 
width, although the KH-unstable ripples are missing. The 
bubbles are also somewhat less collimated, as the material 
close to the Galactic plane is able to cool down more effi- 
ciently and collapse to higher densities. In general, however, 
the bubbles look rather similar. This shows that our results 
are somewhat insensitive to the details of gas heating and 
cooling. This findin g is also cons istent with the an alytical 
predictions of iKina (|2010aj) and King et al.l (120111 ) , where 
by construction an energy-driven outflow occurs when cool- 
ing of the shocked wind and the shocked ISM becomes in- 
efficient; we comment on this result further in Section 15.31 
The effect of the X-ray heating from Sgr A* on the CMZ 
material is even smaller than that on the diffuse ambient 
gas, as explained above. 
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We point out that this insensitivity to quasar photo- 
ionisation heating and radiative cooling is to be expected 
due to the low density and a rather short duration of 
Sgr A* outburst compared with typical cosmological con- 
ditions IjSazonov et al.l 120051 ). At higher densities we would 
expect the structure of the ambient gas to be significantly 
dependent on the details of the cooling function employed 
in the simulations. 



5 DISCUSSION 

5.1 Summary of simulation results 

In general, our numerical simu l ation s appear to confirm 
the suggestion of lZubovas et al.l (|201ll ) that an Eddington- 
limited outburst of Sgr A* outflow is a promising way 
of explaining the morphology of the observed Fermi bub- 
bles. Starting with the model for AGN feedback devel- 
oped to explain the Mbh — Q" relation for classical bulges 
and elliptical galaxies l|King|2003l . 120051 ). the "typical" ini- 
tial conditions for numeri cal simulations of AGN feedback 
(jNavakshin fc Powerll2010r i needed to be ammended only to 
account for (a) the lower present day gas content of the 
Milky Way (i.e., smaller "gas fraction" / g ); (b) the presence 
of a massive gas disc, the CMZ, in the plane of the Galaxy 
in the central 200 pc; (c) a finite duration of Sgr A* out- 
burst, which is a free parameter of our model, and is small 
compared with a dynamical time in a host galaxy bulge. 

Given this setup, we confirm that a spherically symmet- 
ric outflow from Sgr A* is collimated by a geometrically thin 
CMZ disc in directions perpendicular to the Galactic plane. 
The outflow then produces two teardrop shaped cavities that 
have sizes similar to the observed 7-ray emission features. 
We varied the free parameters of our model - the quasar 
outburst duration i q and the gas fraction / g - to constrain 



their values to tq « 1 Myr and / g 



10- 



The former is 



plausible and has interesting implications for Sgr A* feeding 
(cf. i]5.7l below). The latter is poorly c onstrained obs erva- 
tional f y but consistent with estimates bv lMcKed |l990j) and 
ISofud (|201ll ). 

The opening angles of the bubbles as seen from Sgr A* 
are Q,b ~ 0.4 x 4tt. This is significantly smaller than the 
solid angle not obscured by the CMZ in the 'Base' simula- 
tion (fi cmz ~ 0.8 x 47r), showing that CMZ casts a larger 
"feedback shadow" than could be expected based on its ge- 
ometrical aspect ratio alone. This is driven by the following 
two effects. First of all, the CMZ not only hinders the out- 
flow propagation directly through it, but also "reflects" part 
of the outflow. The thermally driven outflow of hot gas ab- 
lated from the CMZ surfaces away from the Galactic plane 
redirects the gas fiowlines towards vertical directions. Sec- 
ondly, the bubbles rise due to buoyancy, and cooler material 
streams to fill the void along the Galactic plane, further re- 
ducing the opening angle of the bubbles. 

In addition, although we were not originally interested 
in the evolution of the CMZ due to Sgr A* feedback - the 
role of the CMZ in our simulations was to stop and redi- 
rect the quasar wind only - we found several CMZ-related 
results interesting from an observational point of view: (a) 
re-shaping of the CMZ into a ring-like structure, perhaps 
explaining Herschel observations that CMZ is a dusty ring 



rather than a disc; (b) an induced star formation mode in 
the CMZ which is trigerred by Sgr A* feedback. These points 
are further discussed in JJ5.6I 

We now make a detailed comparison of our results to 
observations. 



5.2 Gas mass within the bubbles 

In our simulations, two cavities corresponding to the ob- 
served Fermi bubbles are filled with hot and diffuse gas. 
The observationally es timated mass of gas within the bub- 
bles (M bub ~ 10 8 M m ; ISu et al"]|2010l ) is much greater than 
the mass of the shocked wind in our model. The latter, from 
analytical arguments in Section [2] is M w ~ MEddtq ~8x 
10 4 Mo for t q = 1 Myr. The mass contained inside the cavi- 
ties is Mbub, sim ~ 6 x 10 Mq for the 'Base' simulation. The 
density of gas within the bubbles (nbub ~3x 10~ 5 cm -3 ) is 
also much lower than the one typically ado pted in spectral 
modelling of the bubble em ission (n ~ 10~ 2 . ISu et alj|20ld : 
ICrocker fc Aharonianll201ll ). 

We believe that our simulations underpredict the gas 
density inside the bubbles which would otherwise be ob- 
tained in a more sophisticated simulation. We use a one- 
phase model for the ambient medium, whereas observa- 
tions of gas in the Milky Way show that it is multi-phase 
IjDame et al.ll200ll ). We would thus expect some cold and 
warm gas to be present even at large heights above the disc 
before Sgr A* "turns on". Much like in supernova shocks 
expanding into the ambient medium, molecular clouds em- 
bedded in the ambient gas are expected to be overtaken 
by the quasar wind and later evaporate inside the bub- 
bles, increasing the ho t gas density there at late times (e.g. 
iMcKee fc CowieHl975T ). In addition to that, having strong 
density inhomogeneities in the ambient m edium should pro- 
voke strong Ray leigh- Taylor instabilities (King 2010bf) dur- 
ing the quasar outburst event, which we do not model here. 
These may allow formation of dense filaments resilient to 
Sgr A* feedback that are left behind the shock just like the 
molecular clouds discussed above. We hope to explore these 
ideas with improved simulations in the future. 



5.3 Bubble energy content 

The total energy input by the quasar into the system is 



E- m 



I L Ed d t q ~ 8 x 10 56 t 6 erg, 



(14) 



where te is quasar activity time in Myr and r//2 — 0.05 is 
the coupling efficiency of the energy-driven wind (eq. |4|. As 
mentioned in fj2j Ei n is significantly larger than the observa- 
tional estimate of the total gas energy content of the bubble: 
£ bu b,obs ~ 10 54 " 55 erg (|Suet al.ll2010l) . However, in our sim- 
ulations, the injected energy E- m is split between the bubble 
and the surrounding shell. Even when the quasar is active, 
analytical calculations show th at only 1/3 o f the input en- 
ergy remains inside the bubble (jZubovas fc K ing 2012]), and 
this value should decrease after the quasar has switched off, 
as the bubble cools due to expansion. 

To quantify this further, we computed the fraction of 
energy injected into the ambient gas by Sgr A* and retained 
by the bubble in the 'Base' simulation as a function of time 
(other simulations show similar results). To this end, we first 
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Figure 10. Time evolution of the fraction of energy input into the 
system by the quasar wind that is retained inside the bubbles. The 
rest of the energy is contained in the shell of shocked expanding 
ambient gas. 



compute the sum of kinetic and thermal energy for all the 
gas particles in the simulation, E to t(t). We exclude in this 
calculation the gravitational potential energy of gas because 
for most SPH particles in the calculation changes in the 
latter are small compared with the changes in the kinetic 
and thermal energy. We then define the energy content of 
the gas within the bubble, -Ebub(t), as the sum of kinetic and 
thermal energy of the gas inside the bubble. The fraction of 
energy retained in the bubble is thus defined as 



£B,in(t) = 



A-Ebub _ Ebub(t) — £'bub(0) 

AEto, ~ 



Et, 



(15) 



t(t)-E t ot(0) 

Figure [TTJ presents the evolution of 6E,in with time. 
While the quasar is active, the fraction of energy retained 
inside the bubble grows slightly and fluctuates around the 
analytically derived value of 1/3. Once the quasar switches 
off, the expanding bubble cools and transfers most of its en- 
ergy to the surrounding shell, so that by t = 6 Myr, only 
~ 3% of the total energy input is retained inside the bubble. 
Hence the total energy content of the bubbles by the end of 
the simulation is 



E h 



0.03£ in ^ 2.5 x 10 55 t 6 erg, 



(16) 



a value similar to but still a little larger than the observed 
one. 

A further decrease in the bubble energy content may be 
caused by radiative cooling. To check this, we consider the 
evolution of E to t (£) and find that after t — i q it is conserved 
to within a few %. Adding a reasonably parametrized (us- 
ing r s — 1 kpc as a scale radius) value for the gas potential 
energy does not change the result either. Therefore the im- 
portance of cooling on the total energy content is minimal 
and does not affect our previous considerations. 



5.4 Expected radiation from the bubble 

The main emission compo nents of intere st are the 7-ray 
emission from the the lobes (|Su et al.ll2010l ). a tentative mi- 
crowave feature coinci dent with the bubbles identified in the 
WMAP all-sky maps l|Finkbeinerll2004l ; ISu et aIll201Ch and 



an X-shaped feature visible in the X-rays closer to the Galac- 
tic plane, coinciding wit h the edges of the bu bbles, known 
as the "ROSAT limbs" (|Snowden et ail 1 19971 ) . We discuss 
each of the three components in turn. 

While we do not model cosmic ray (CR) particles in our 
paper, we note that astrophysical shocks are known to accel- 
erate electrons, protons an d other particles to CR energies 
iJBlandford fe E ichlcr 1987) in a variety of environemnts and 
put as much as ~ 10% of the blast wave energy into the high 
energy particle component in the case of supernova shocks. 
We certainly do not see an obvious reason why shocks driven 
by an even faster outflow from a quasar would be less effi- 
cient in producing CRs than supernovae. 

Two competing explanations for the origin of the 7- 
ray emission from the Fermi bubbles wer e sugge sted in the 
literature to date. ICrocker fc Aharonianl (|2011T ) suggested 
that the emission is powered by CR protons through pp 
collisions with the plasma in the bubbles. In this scenario, 
~ 10 39 erg s _1 of energy is injected into the bubbles in the 
form of CRs for ~ 10 10 yrs to achieve saturation in the sys- 
tem (cf. also ICrockeil 1201 ll i This yields ~ 3 x 10 58 erg in 
CR enegy alone. This is some 3 orders of magnitude larger 
than the thermal energy retained within the bubbles in our 
"Base" simulation fequation ll6p . Therefore we conclude that 
our model is very unlikely to produce a bright enough 7-ray 
emission if emitting parti cles are hadrons. 

On the other hand, iMertsch fc Sarkan (|201lh explore 
the possibility that electrons accelerated by shocks in tur- 
bulent plasma inside the bubbles can reproduce the spectral 
features, as well as the constant surface brightness profile. 
They find that this acceleration process can continuously 
replenish the energetic electron population, overcoming the 
problem of rapid electron cooling via the inverse-Compton 
process. The cooling time due to both IC and synchrotron 
losses is less than 5 Myr for 100 GeV electrons at z = 5 kpc 
abov e the Galactic plane (Fig. 28 and Section 7.1 in lSu et ail 
120101 ). Since the current 7-ray luminosity of the bubbles is 
~ 4 x 10 37 erg s _1 , and the age of the bubbles is 6 Myrs in 
our model, we conclude that at least ~ 10 52 erg of CR elec- 
trons is required to explain the Fermi Bubbles, and some- 
what more if we account for continuous replenishment from 
plasma turbulence. This would be a very small fraction of 
the bubble energy content in our model, and on that basis 
at least is possible. 

The origin of "continuous" shock driving inside the bub- 
bles can be naturally explained by our model. Since the 
Fermi bubbles are still expanding at the end of the sim- 
ulations, we expect that the whole volume of the cavities 
inflated by the quasar outburst is still filled with shocks. 
As noted in JJ5.2I a higher resolution modelling that in- 
cludes a multi-phase description for the ambient pre-shock 
gas is likely to result in strong Rayleigh- Taylor instabilities 
of the wind shock front during the quasar outburst event 
l|Kinel l2010bl i This should let cold gas filaments fall into 
the bubble's interior due to gravity, leading to high Mach 
number shocks. Protons and electrons could thus be con- 
tinuosly accelerated on these shocks now, as suggested by 
IMertsch fc Sarkarl (|201ll 1. 

The observed microwave haze spatially coincides with 
the Fermi bubbles, although its intensity decreases with 
height above the Galactic plane, especially in the southern 
bubble (see the bottom right panel of Figure 18 in lSu et ail 
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|201y) . The luminosity of synchrotron emission from cooling 
electrons decreases with height above the plane possibly due 
to an expected decrease in the magnetic field strength. We 
believe that same electrons responsible for the 7-ray emis- 
sion may be responsible for the WMAP haze. 

The ROSAT-observed X-ray background in the region 
of the Fermi bubbles is composed of X-shaped ridges coinci- 
dent with the lower limbs of the bubbles, most clearly visible 
in the hard (1.5 keV) band of the instrument, and a cavity 
in soft-band X-ray emission coinc i dent with the bulk of the 
bubble area l|Snowden et al.ll 19971 ; ISu et al.ll2010T ). The spa- 
tial coincidence of the features suggests a common origin. 

In our simulations, gas inside the bubbles has temper- 
atures as high as ~ 10 — 20 keV (see Figure [H right) and 
so most of its thermal emission is harder than what ROSAT 
would detect. This could then explain the deficit of soft X- 
ray emission from within the bubbles. The shocked ISM (see 
Section \5 . 5 1 below) is heated to temperatures of a few times 
10 7 K, (although that is expected to scale as ecu 2 , so may 
vary somewhat if the model for the potential in the inner 
10 kpc of the Galaxy is varied). The shocked shells around 
Fermi bubbles should then be most visible in the harder 
ROSAT bands. X-ray emissivity scales as oc n 2 , where n is 
gas particle density. The latter tends to be larger close to 
the Galactic plane due to superposition of the shells around 
the two bubbles (Figure [3] right) and vertical stratification. 
Therefore we may expect the lower limbs of the shells to be 
brighter in X-rays, as observed. 



5.5 Bubble edges and the outer ISM shock 

The expanding SMBH wind drives a strong forward shock 
into the ISM, whic h is expected to be adi abatic and move 
with wism ~ 4/3t>e (|Zubovas &: King|2012i ). Since the shell 
thickness is initially zero, its vertical extent should also 
be nsM ~ 4/3rb, i.e. the thickness of the snowplough 
shell should be l/3rd of the bubble height (jZubovas fc Kind 
|2012P |. We note that in the simulations, it is somewhat 
greater: d s h/?"b ~ 0.4. This result is not due to poor resolu- 
tion. At late times, the thickness of the shell is ~ 10 times 
greater than the typical SPH particle smoothing length in- 
side it (ftsPH.sh ~ 500 pc), therefore we believe the large- 
scale morphology of the region is well resolved. 

The most likely reason for the discrepancy is the fact 
that the analytical model assumes a continuous quasar out- 
flow driving the shell, whereas in our simulations the quasar 
turns off at t — 1 Myr. As a result the shell is less compressed 
in the radial direction. 

In addition, our one-phase treatment of the ambient 
diffuse medium (cf. i]5.2l) can under-estimate the radiative 
cooling in the shell. Therefore, a more realistic multi-phase 
simulation could be expected to have thinner shell envelop- 
ing the Fermi bubbles. 

Simulations also show transition regions at the edges of 
the bubbles ( e.g. at r cy i = 4—5 kpc in Figure(4] right). These 
are < 1 kpc thick zones where both temperature and density 
change from the values appropriate for bubble interior to 
those of the swept-up ISM. The observed bubble edges are 
5 — 10° wide, corresponding to a few hundred parsecs, in 
good agreement with simulation results. However, we must 
note that the SPH particle smoothing length in these regions 



is comparable to the region thickness, therefore numerical 
effects probably dominate the result. 



5.6 Has Sgr A* feedback affected the CMZ? 

The primary goal of our paper is investigating whether 
Sgr A* feedback is a reasonable model for the morphology 
of the Fermi bubbles. In doing so we purposely introduced 
a very simple model for the CMZ - a £(-R) oc R~ x cir- 
cular disc in the plane of t he Galaxy with mass an d size 
consistent with observations (jMorris fc Serabvnlll996r ). This 
simple model nevertheless resulted in an interesting trans- 
formation and exitation of the CMZ that is worthy of further 
discussion. 



5.6.1 The Herschel ring: a feedback-compressed disc? 

In particular, we find that, under the strong coercion from 
Sgr A* feedback, the initial disky configuration of the CMZ 
attains a morphology more reminiscent of a dense ring (see 
Section 14.1.31 and Figure 0). We believe this is a general 
result of an AGN feedback acting on a disk since the physics 
behind this "anti-diffusion" evolution of the disk is simple 
and thus robust. 

When the inner regions of the disc are blasted with 
an outward-directed feedback from the centre, little of the 
material outflows to infinity; most actually falls back onto 
the disc at large radii, mixes with gas there and thus induces 
an inward-directed radial flow at those radii. AGN feedback 
thus serves as an external agent that forces shock mixing of 
gas initially located at different radii and carrying different 
specific angular momenta. An initially broad distribution of 
angular momentum (a disc) becomes more narrow (a ring). 

The slow viscous evolution of the dense ring that is 
formed by feedback suggests that such a ring should be still 
present in the GC today. In fact, recent Herschel observa- 
tions of the GC region have revealed a ring -like structure 
on the scale of ~ 100 pc (|Molinari et alj|2011f ). qualitatively 
similar to what we see in the simulations. The observed fea- 
ture is elliptical and offset from Sgr A* . In Section 15.6.21 
below we discuss how this strongly non-circular and offset 
structure could be formed in our model. 

The fact that the CMZ is only perturbed but not dis- 
persed to infinity is consistent with analytical predictions in 
H2.3I We note also that it may seem paradoxical how this 
very massive gas disc gets evacuated in the inner region and 
yet Sgr A* continues to accrete from presumably a much 
smaller disc at the assumed Eddington rate for ~ 1 Myr as in 
our "Base" simulation. However, the paradox is easily solved 
by realising that weight per unit mass of the gas - i.e., accel- 
eration due to gravity - scales as oc M(R)/R 2 , where M(R) 
is the total mass enclosed within radius R. The disc that fed 
Sgr A* may have ha d a radius of orde r R ~ 0.1 pc or even 
less JNavakshin fc Cuadral2005l ; I Alexander et £01 20121. The 
weight of the accretion disc of mass M^ = 1O 5 MsM0 is then 

Wac = GMbhAfd = L2 x 10 36 Ms dyn ( (17) 

R 
which is some ~ 20 times larger than that estimated for 
the CMZ (equation [5]). It is thus perfectly reasonable that a 
smaller-scale disk would be too tightly bound to the SMBH 
to be expelled or even significantly perturbed by feedback, 
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unlike gas at larger radii. Nayakshin, Power & King (2012, 
Ap J submitted) show in their Appendix that this is a general 
point. 



5.6.2 Non-circular orbit of the Herschel ring 

iMolinari et al.l (|201ll ) find that the Herschel ring is offset 
from Sgr A* and thus the centre of the Milky Way by about 
50 pc, which is comparable with the size of the ring itself (see 
their figure 5). Any offset, and especially such a large one, is 
not naturally expected if the ring is a long-lived structure. 
While our model shows no offset, we note that there may be 
a natural way to produce that, although further numerical 
simulations are needed to confirm these ideas. 

In particular, we have assumed here that the feedback 
from Sgr A* is exactly spherically symmetric. This is likely 
to be an over-simplification. Even if wind outflow is quasi- 
spherical on small scales (e.g., hundreds of Schwarzchild 
radii) , on the somewhat larger scales of stellar discs we know 
that there must be a an accretion disc which cannot be read- 
ily expelled by Sgr A* feedback as discussed above. We also 
know that the young stellar discs orbiting Sgr A* are not 
simple planar structures, with signific ant warps needed t o 
explain the observed stellar kinematics (jBartko et al.l l2009). 
Large warps also natural ly occur in nume rical simulations of 
star forming gas flows (Hobbs & Navakshin 2009). There- 
fore, due to shadowing of Sgr A* feedback by the ~ 0.1 
pc gas flow (which may be strongly non-circular), one side 
of the CMZ may have experienced a different amount of 
feedback compared with the other side, causing a strong 
non-axysimmetric perturbation to the CMZ disc and per- 
haps leading to an offset ring reminiscent of the observed 
Herschel ring. We note that the asymmetry in the feedback 
should not be too large over the whole 47r solid angle, how- 
ever, so as to not result in too asymmetric Fermi Bubbles. 

Additionally, natural non-axisymmetric density varia- 
tions in the pre-feedback CMZ may result in outflow within 
the CMZ proceeding at different velocities in different direc- 
tions, producing the offset. 



5.6.3 Induced star formation in the CMZ? 

A particularly interesting outcome of the CMZ compression 
by the outflow is the formation of a dense clump of gas at t ~ 
3.5 Myr (see Section l4T3|) . The large mass, ~ 10 7 M Q , of the 
clump and its position at R ~ 100 pc from the Galactic cen- 
tre sugg est that the observed young Galactic ce ntre stellar 
clusters l|Serabvn et al.lll998l ; ifiger et al.lll999allbh or some 
of th e largest molecula r clouds (e.g. Sgr B2, IScoville et al.l 
1 19751 ; iReid et al.l 12009m may have formed this way. Many 
smaller clumps in the dense ring may be susceptible to star 
formation; recent I R observations (|Yusef-Zadeh et al.ll2009i ; 
llmmer et al.l |2012j) reveal that the star formation rate in 
the CMZ during the past ~ 1 Myr was ~ O.O8M yr" 1 . 
At this rate, all of the CMZ gas would be turned into stars 
in ~ 6 x 10 8 yr; even if the CMZ was supplied by molecu- 
lar gas from further out, the whole molecular gas content 
of the Milky Way would have been turned into stars in 
~ 6 Gyr. This suggests that the current star formation rate 
in the CMZ is higher than the long-term average, further 
implying a rather recent perturbation, consistent with the 



results of our simulations. Furthermore, combined Chandra 
and HST observations have revealed a population of isolated 
massive (O and B gian t, but also WR) stars in the CMZ (e.g. 
iMauerhan et alj|2010h . These stars are not associated with 
any of the known clusters and may have formed in small 
associations. Our simulations suggest that while these stars 
are isolated now, they may have formed coherently in time 
if not in space due to a single very powerful perturbation 
produced by Sgr A* quasar outburst that turned the disc 
into a star-forming ring. 

Finally, we note that the strong perturbations to the ini- 
tially circular orbits of gas in the CMZ by Sgr A* feedback 
could form gas clumps on eccentric orbit. Indeed, the orbit of 
the clump found in the "Base" simulations is mildly eccen- 
tric, e ~ 0.2. That was obtained in a perfectly azimuthally 
symmetric simulation; any deviation from this assumption 
would have likely resulted in an even more eccentric clump. 
Such eccentric clump formation mechanism may be relevant 
to the origin of the Arches cluster which has a rather non- 
circular orbit (|Stolte et al.ll2008l ). 



5.7 Implications for AGN feeding models 

Sgr A* is the closest SMBH, and while its dimness is 
a familiar (and important) tale that stimulated develop- 
ment of non-radiative models of ac cretion for low den- 
sity gas flows near SMBH s (e.g., iNaravan et al.l Il995l ; 
iBlandford &: Begelmanlll999i ). little has been known about 
the past of Sgr A* as an accretion-powered SMBH. This 
is mainly due to potential difficulties of discovering signs 
of past Sgr A* activity. Light-echos due to X-ray fluores- 
cence of Sgr A* radiation on molecular clouds is a power- 
ful technique (|Sunvaev fc Churazovlll99Sl ; iRevnivtsev et al.l 
12004 ; iTerrier et alj|20ich but can only be used to constrain 
Sgr A* activity up to ~ 10 4 yrs at best due to the finite size 
of the Galaxy and the ~ 1/R 2 fading of signals from clouds 
at large distance R from the Galactic Centre. 

However, shocks induced by an outflow from Sgr A* 
may continue to be visible for about a dynamical time of 
the Galaxy, i.e., Rg/o ~ 50 Myrs, where we set Galaxy 
"radius" to be ~ 5 kpc and velocity dispersion a — 100 
km s _1 for il lustrative purpo ses. Therefore, the Fermi lobes 
detected by ISu et al.l i|2010f ) may be such a "shock-echo" 
of the past Sgr A* activity. Following IZubovas et al.l (J201l|) 
we required the activity episode of Sgr A* to coincide with 
the well known star f ormation event in the Galactic Centre 
~ 6 million years ago IjPaumard et al.ll2006T ) , and found that 
the model has a number of attractive observational conse- 
quences which lend some support to this picture. 

If Sgr A* indeed did have an Eddington-limited out- 
burst for as long as ~ 1 Myr, it must have accreted ~ 1O 5 M0 
of gas. This is an o rder of magnitude more than the mass 
in th e stellar discs l|Paumard et al.l 120061 ; iNavakshin et al.l 
I2006J ). The fact that Sgr A* managed to accrete ~ 90% of 
the gas from the star-forming accretion disc is highly sig- 
nificant for the general question of how SMBH are fed, and 
should be explored further as a potential example of an AGN 
disc that avoided "star formation catasrophe" in which gas 
is believed to be turned into stars too r apidly to feed AGN 
|Goodmanll2003l ; INavakshin et alj|2007h . 
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5.8 Uncertainties and deficiencies of our work 



As discussed in £15.41 we model the ambient gas "halo" with 
a single phas e medium, w hereas the ISM is expected to be 
multi-phase (|McKeeHl990T ). The lower density medium, es- 
sentially unresolved in our simulations, is probably quite 
important for the dynamics of the fast outflow from Sgr A* 
interacting with the gas in the bulge of the Galaxy (see 
also i]5.2l) . Therefore some of our conclusions (e.g., the likely 
quasar phase duration t q , the geometrical thickness of the 
shocked shell, etc.) may somewhat depend on the treatment 
of the ambient medium, and future work is needed to quan- 
tify this issue. 

Also, a more realistic model for the geometrical arrange- 
ment of the diffuse "halo" gas that the outflow interacts with 
is desirable as this probably affects the eventual shape of the 
bubbles. In our simulations, we consider a spherical halo 
mass distribution, whereas in reality, we expect stratifica- 
tion on a large scale, with higher gas density in the Galactic 
plane than perpendicular to it, and higher density in the 
bulge than outside. We can make a very simple estimate of 
the magnitude of this effect by considering the dependence 
of outflow stalling radius on the gas fraction, which we take 
here to be the mean gas fraction along a direction of ex- 
pansion. The analysis in Paper I and Section [2] shows that 
■Rstaii oc v% oc / g . As the bubble expands in an approxi- 



mately self-similar fashion (see Section [47TJ) after the quasar 
switches off, we can estimate that in any given direction 
Rbub oc / cft . Therefore a factor of three difference in gas 
fraction results in approximately a factor of two difference in 
bubble radius. We hope that future observations of diffuse 
gas in the bulge of the Milky Way will help to constrain this 
part of our model. 

Our results are, as expected, sensitive to the total 
mass of the CMZ, producing bubble morphology inconsis- 
tent with observations if it is reduced to 1O 7 M0 or less. 
This low value of the CMZ mass is however unrealistically 
small, smaller than the current observational constraints o f 
3-5xl0 7 M© (JDahmen et al.lll998l : lpTerce-Price et"al1l2000h . 
A potentially important caveat is that we model the CMZ 
as a smooth disc, rather than as a clumpy distr ibution of 
gas iJMolinari et all l201ll : [Morris fc Serabvnlll996l '). In prin- 
ciple, a real outflow may stream past the clumps and so 
the collimation effect of the CMZ would be smaller than 
what we find. However, observations also show that the col- 
umn density of gas toward the Ga lactic centre does not vary 
strongly with viewing direction IjGoto et al.ll2008l ). so the 
CMZ should provide a strong covering effect to any outflow 
coming from Sgr A*. In addition, we find that the CMZ as- 
pect ratio has a negligible effect on the bubble properties, 
i.e. there is a range of CMZ column densities for which the 
outflow is collimated efficiently. Therefore, we believe that 
the uncertainty in our results due to the shape of the CMZ 
is relatively small. 



6 CONCLUSIONS 

We have presented numerical simulations of a wide angle 
outflow from Sgr A*, temporally coincident with the star 
formation event ~ 6 Myr ago, and collimated into directions 
perpendicular to the plane of the Galaxy by the presence of 



a massive disc of gas, the Central Molecular Zone. Our main 
results are: 

• This model is a plausible way to inflate the 7-ray emit- 
ting lobes recently observed by the Fermi-LAT, provided 
that Sgr A*'s outburst duration is ^S 1 Myr. 

• The energetics of the model is consistent with the ob- 
served 7-ray emission if radiating particles are electrons 
rather than hadrons. 

• Sgr A* feedback could have reshaped the CMZ from a 
disc-like configurati on into a ring-like on e reminiscent of the 
observed structure IjMolinari et al.ll2011i ) . We speculate that 
somewhat assymetric feedback could produce an observed 
offset of the ring and eccentric orbits for young star clusters 
such as the Arches star cluster. 

• Furthermore, CMZ compression by the feedback out- 
flow may explain the formation of dense GMCs, the young 
star clusters, the population of isolated massive stars, and 
also result in the present-day high star formation rate in the 
Galactic Centre. 

An important side result of this paper is that the same 
feedback model appears to work for quasars as their SMBHs 
establis h the M — a relation and clear the host ga laxies of 
gas fe.g. lKindl2lMl2005l ; lNavakshinfe Powerll2010I ). as well 
as for a short burst of activity in Sgr A*, a SMBH that is 
somewhat below the M — a relation in a quiescent galaxy. 
This finding suggests that there is nothing fundamentally 
different between feedback from SMBHs at gas-rich epochs 
(z > 2) and that from local galactic nuclei, except for the 
lower amount of fuel they receive. 

More detailed future treatment of the feedback process, 
with improvements in both the physics of the simulations 
and more realistic observationally constrained initial con- 
ditions, may provide interesting and unique constraints on 
cosmological models of AGN feedback. 
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